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We present recent developments on the role of the running coupling constant at the 
intersection of perturbative and nonperturbative QCD. A number of experiments show 
a smooth transition from small to large scales given by the four-momentum transfer in 
the reactions. This is at variance with perturbative QCD where the running coupling 
constant becomes infinite when the scale equals Aqcd- Approaches using an effective 
coupling constant could help interpret the opposite trend of data as compared to standard 
perturbative QCD predictions. We give an overview of the role of the coupling constant 
in the procedure to match nonperturbative hadronic model to perturbative QCD and 
we propose an extraction of an effective coupling constant from inclusive electron proton 
scattering data at large Bjorken x. 
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1. Hadron Phenomenology 

The description of the internal structure of the strongly interacting particles is one 
of the key goals of QCD. At moderate energy scales, the hadronic representation 
supersedes the partonic description, making it challenging to describe the dynamics 
of scattering processes and hadronic structure. On the other hand, it is well known 
that the Parton Distribution Functions (PDFs) present a framework for connecting 
the low and high-energy regimes. Deep Inelastic processes allow us to look with a 
good resolution inside the hadron and to resolve the very short distances, i.e. small 
configurations of quarks and gluons. The insight into the structure of hadrons is 
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regulated by factorization theorems: the large virtuality of the photon, Q 2 , involved 
in DIS processes allows for the factorization of their amplitudes into hard (pertur- 
bative) and soft (nonperturbative) contributions. At short distances, one singles out 
a hard scattering process described through Perturbative QCD (PQCD). The large 
distance part of the process, i.e. the PDF, reflects how the quarks and gluons are 
distributed inside the target. 

So far, we have introduced two different concepts: a first transition from the 
hadronic to partonic representation is defined by a hadronic scale of a few hundred 
MeV 2 (see e.g. discussion in Ref. [I]); a second scale, namely, the factorization scale, 
typically > 1 GeV 2 , is introduced through the virtuality of the photon. 

Although, as we have just explained, the perturbative stage of a hard collision 
is distinct from the nonperturbative regime characterizing the hadron structure, 
experimental observations suggest that, in specific kinematical regimes, both the 
perturbative and nonperturbative stages arise almost ubiquitously, in the sense that 
the nonperturbative description follows the perturbative one. This third concept is 
known as parton-hadron duality, and we will here understand it as being yet another 
manifestation of the perturbative to nonperturbative transition in QCD. 

1.1. Can we extract a s in the infrared regime from hadronic 
phenomenology ? 

The standard procedure to fix the hadronic (nonperturbative) scale pushes pertur- 
bative QCD to its limit. The hadronic scale turns out to be of a few hundred MeV 2 , 
where the strong coupling constant has already started approaching its Landau pole. 
However, the relative stability of the N m LO evolution is what justifies the perturba- 
tive approach. An even lower hadronic scale makes sense when considering a freezing 
of the coupling constant in the infrared region. The standard scheme of scale fixing 
can then be extended to a nonperturbative evolution framework where the effective 
coupling is free of Landau pole, e.g. Refs. [2H1], while being parameterized with a 
physical set of parameters [6j [7] . 

On the other hand, the factorization scale fixing procedure is carried out entirely 
within the domain of PQCD, thus relying on the knowledge of a s . Although these 
two approaches have been considered so far complementary to each other, a unified 
description might derive through the definition of the effective coupling, as they 
both broaden the ways of analyzing the freezing of the running coupling constant. 
It is in this direction that the new procedure proposed in Refs. [6l [7] broadens 
the ways of analyzing the freezing of the running coupling constant: T-odd TMDs 
are possible candidates to study the behavior of a s at intermediate and low Q 2 , 
along with the observables proposed in Ref. [8] (polarized DIS sum rules) and in [9] 
(perturbative evolution of large x proton structure functions) . 

In this contribution we focus on the analysis of Ref. [9] , where the implications 
of parton-hadron duality are explored in the large x region of inclusive electron 
proton scattering experiments (Bloom-Gilman duality |10j). Our ultimate goal is 
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to provide a procedure on how to match nonperturbative models to PQCD, using 
experimental data. The relevant kinematical variables are: x = Q 2 jIMv (M being 
the proton mass and v the energy transfer in the lab system), the four-momentum 
transfer, Q 2 , and the invariant mass for the proton, P, and virtual photon, q, system, 
W 2 = (P + q) 2 (W 2 = Q 2 (l/x - 1) + M 2 ). For large values of Bjorken x > 0.5, 
and Q 2 in the multi-GeV 2 region, one has W 2 < 5 GeV 2 , i.e. the cross section 
is dominated by resonance formation. While it is impossible to reconstruct the 
detailed structure of the proton's resonances, these remarkably follow the PQCD 
predictions when averaged over x. Although Bloom-Gilman duality was observed at 
the inception of QCD, quantitative analyses could be attempted only more recently, 
having at disposal the extensive, high precision data from Jefferson Lab [TT]. PQCD- 
based studies [9l IT2UT4] . have been presented that include higher- twist contributions 
or, more generally, the evidence for nonperturbative inserts, which are required to 
achieve a fully quantitative fit. In addition, NLO PQCD evolution at large x can 
be sensitive to Large x Resummation (LxR) effects. The consequence of LxR is a 
shift of the scale at which a s is calculated to lower values, with increasing x (see 
for instance Refs. [HI US]). This introduces a model dependence within the PQCD 
approach in that the value of the QCD running coupling in the infrared region being 
regulated by LxR (as we explain below) simultaneously leads to a suppression of 
higher- twist effects. The higher- twist effects get, in fact, absorbed in the coupling's 
infrared behavior. 

In Section 2 we single out the nonperturbative outcome from a quantitative anal- 
ysis of the Bloom-Gilman parton-hadron duality. Specifically, we analyze the role 
of the running coupling constant in the infrared region in tuning the experimental 
data Q21 [T7] . The new approach on the freezing of the running coupling constant 
and its role on the extraction of TMDs from experiment [51 [7] is summarized in 
Section 3. In Section 4 we draw our conclusions. 



2. From perturbative to nonpertubative QCD 

Bloom-Gilman duality implies a one-to-one correspondence between the behavior 
of the structure function, F%, for unpolarized electron proton scattering in the reso- 
nance region, and in the PQCD regulated scaling region. A quantitative definition of 
duality is accomplished by comparing limited intervals (integrated in Bjorken- a; over 
the entire resonance region) defined according to the experimental data, e.g. [HI 
119] . Hence, we compare the scaling results as a theoretical counterpart, or an output 
of PQCD, in the same kinematical intervals, and at the same scale Q 2 . Namely we 
consider the ratio, 

J Xmln dxF» ls (x,Q 2 ) 

We consider duality to be fulfilled if i? Res / DIS (x avc , Q 2 ) is 1. Note that the defini- 
tion in Eq.([T]) relies on the fact that the PQCD evaluation, F^ IS (a;, Q 2 ), is very well 
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constrained in the region of interest (x > 0.3) despite it does not correspond di- 
rectly to measured data. F^ is an input that once fed into the evolution equations 
determines the structure functions behavior at much larger Q 2 . However, the error 
on this type of backward evolution is expected here to be small, being dominated 
by the valence contribution (a quantitative analysis of the latter will be carried out 
in an upcoming study [H]). Had we applied the same procedure to low x where the 
singlet and gluon distributions govern F2, we would have gotten a much larger error 
at low Q 2 because of the strong correlation with the value of a s . 

Besides perturbative evolution one has to take into account several hadronic 
corrections to -F^ IS . For instance, if we evolve the structure functions to NLO, we 
find that duality is violated by a given amount. However Target Mass Corrections 
(TMCs) are important here and move the ratio closer to unity. The most important 
effect for our purposes is the effect of LxR, that we develop hereafter. 

2.1. Large-x Resummation 

Large x threshold resummation effects (LxR) arise formally from terms containing 
powers of ln(l — z), z being the longitudinal variable in the evolution equations, that 
are present in the Wilson coefficient functions C(z). Below we write schematically 
how the latter relate the parton distributions to e.g. the structure function F2, 



where we have considered only the non-singlet (NS) contribution to F2 since only 
valence quarks distributions are relevant in our kinematics. The logarithmic terms 
in C(z) become very large at large x, and they need to be resummed to all orders 
in a s . Resummation was first introduced by linking this issue to the definition of 
the correct kinematical variable that determines the phase space for the radiation 
of gluons at large x. This was found to be W 2 = Q 2 (l - z)/z, instead of Q 2 [I51I20] . 
As a result, the argument of the strong coupling constant becomes z-dependent: 
a s{Q 2 ) — > cts(Q 2 (l — z )l z ) EU [22]. In this procedure, however, an ambiguity is 
introduced, related to the need of continuing the value of a s for low values of its 
argument, i.e. for z — ► 1 [23] . 

Since the size of this ambiguity is of the same order as the higher-twist correc- 
tions, it has been considered, in previous work [24], as a source of theoretical error 
or higher order effects. We propose an accurate analysis |17j from which one can 
extract a s for values of the scale in the infrared region. To do so, we investigate the 
effect of varying the form of the running coupling on the evolution equations. We 
consider the following choices: 





• a s (Q 2 ) ; 

• an expansion of a s (W 2 ) in ln((l — z)/z), to NLO, 
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Fig. 1. Running of the strong coupling constant in the MS with A L o = 174MeV. The solid black 
curve represent the LO ct s (Q 2 )/47T. The dashed red curve represents the expansion of the strong 
coupling in ln((l — z)/z), for z = 0.7 ; the dotted blue curve is the complete a s (Q 2 (l — z)/z)/4m 
for the same value of z. 



• the complete z dependence of a s (W 2 ). 

The running of each of the three versions of the coupling constant starts being 
very different when z — y 1. We illustrate this behavior on Fig. Q]for z = 0.7. The 
infrared behavior of the coupling constant with argument Q 2 starts to matter to 
lower Q 2 values than for a coupling constant which argument is W 2 . On the other 
hand, the corresponding Landau poles do not coincide: for small Q 2 values and large- 
z, the argument of a s (W 2 ) differs from the logarithmic terms taken into account in 
the NLO expansion. The asymptotic value differs for a s (W 2 ) and expansion w.r.t. 
a s (Q 2 ), as shown in the inner frame. 

The meaning of LxR becomes very clear from Figjl] It is now understood that 
the only free parameter in testing the realization of duality here, is related to a s . 
By playing with the argument of the running coupling constant, we can tune the 
scaling structure function and extract the low Q 2 behavior that determines duality. 
For instance, by setting a maximum value for z one would prevent the DGLAP 
evolution from including extremely large values of the coupling constant. Moreover, 
this z max could define a criterion of convergence of the expansion w.r.t the complete 
a s (W 2 ), as illustrated in Fig.d 

This exercice has to be repeated for each experimental data point. We observe 
from our analysis that the maximum value of z or, equivalently, the scale in which 
the running of the coupling is stopped changes from one to another data point. 
A rough qualitative parameterization of the realization fo duality would look like 
Fig. [3j where we have use the Cornwall's effective charge resulting from a massive 
gluon propagator [2]. In effect, the dynamical gluon mass generation leads to the 
freezing of the QCD running coupling constant. The nonperturbative generalization 
of a s (Q 2 ) comes, here, in the form 
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Fig. 2. Running of the strong coupling constant in the MS with A L o = 174MeV at fixed Q 2 . The 
solid black curve represent the LO a s (Q 2 )/4-K. The dashed red curve represents the expansion of 
the strong coupling in ln((l — z)/z) ; the dotted blue curve is the complete a a (Q 2 (l — z)/z)/4:n. 
The cut of the 



(4) 

The zero gluon mass limit leads to the LO perturbative coupling constant momen- 
tum dependence. The m 2 (Q 2 ) in the argument of the logarithm 

-1-7 

(5) 

with (7) = 1/11, tames the Landau pole. The mass m 2 (Q 2 ) can be understood as a 
constituent gluon mass, and depends non-trivially on the momentum transfer Q 2 . 
As a consequence a s (Q 2 ) freezes at a finite value in the IR [2 [3]. 

An analysis of the complete data set collected at JLab in addition to the existing 
large x data will allow us to carry our study to a quantitative level [17] . Three 
nonperturbative representations of the coupling constant will there be analyzed (2j- 
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Fig. 3. The running coupling constant. The solid red curve represents the perturbative coupling 
constant in MS with A L o = 174 MeV. The blue curve represent the Cornwall effective charge of 
Eq. l|4} for the same value of A, = 1.05A 2 and for p = 1.7, 1.66 on, respectively, the solid and 
dashed curves. 
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3. The other way: From nonperturbative to perturbative QCD 

Now, trying to match nonperturbative model with PQCD evolution, we examine 
the consequences on the hadronic scale in considering the same Cornwall's effective 
charge in the evolution of the PDFs, evaluated in hadronic models. 

In QCD all matrix elements must have a scale associated to them as a result of 
the RGE of the theory. A fundamental step in the development of the use of hadron 
models for the description of properties at high momentum scales was the assertion 
that all calculations done in a model should have a RGE scale associated to it [25] . 
The momentum distribution inside the hadron is only related to the dynamical scale 
and not to the momentum governing the RGE. Thus a model calculation only gives 
a boundary condition for the RG evolution as can be seen for example in the LO 
evolution equation for the moments of the valence quark distribution 



where d r ^ s are the anomalous dimensions of the NS distributions. For calculations 
in the bag model, the dynamics it describes is unaffected by the evolution pro- 
cedure, and the model provides only the expectation value, (q v {^o))n, which is 
associated with the hadronic scale. When considering the nonperturbative solution 
of the Dyson-Schwinger equations, that results in the appearence of an infrared 
cut-off, the gluon mass Eq. ([3]) will only affect the evolution of the PDFs. The gen- 
eralization of the coupling constant results to the structure function imply that the 
LO evolution Eq. (j6)) simply changes by incorporating the nonperturbative coupling 
constant evolution Eq. ((4]). In other words, the hadronic scale that ensues from the 
nonperturbative procedure is quantitatively different from the perturbative scheme. 
However, the results from both procedure are close, what give us confidence on the 
perturbative procedure even at low scales. We note however, that the corresponding 
hadronic scale, for the sets of parameters chosen in Fig. [3] turns out to be slightly 
smaller than in the perturbative case (/Zq ~ 0.1 GeV 2 ). The physical meaning here 
would be that the nonperturbative approach seems to favor a scenario where at the 
hadronic scale we have not only valence quarks but also gluons and sea quarks. 

The application of the nonperturbative framework to the evaluation of the T- 
odd TMDs illustrates the uncertainty on model predictions coming only from the 
matching with RGE. In Ref. {S\ we have considered the errorband resulting from the 
uncertainty of the initial value of a s (/Jo) on the bag model evaluation of the Sivers 
and Boer-Mulders functions. This observation shows that the naive scenario may 
well serve to make predictions, within a reasonably small band, which should not 
be far from experimental expectations. Within this interpretation, T-odd TMDs 
are possible candidates to study the behavior of a s at intermediate and low Q 2 , 
though still biased by strong theoretical difficulties such as the TMD evolution [26] 
or ambiguities due to the choice of the hadronic representation. 
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4. Conclusions 

A combined analysis of the extractions of the running coupling constant in the in- 
frared region suggests a novel definition of the effective charge [17], following the 
example of Ref. [8] where the effective coupling constants are phenomenologically 
inferred from different processes, and to calculations based on Schwingcr Dyson 
equations. This analysis also bears potential important consequences for the con- 
nection of low scale/hadronic models with experiments in the multi-GeV region. 

Acknowledgments 

We are grateful to J. P. Chen, A. Deur, and V. Vento for fruitful discussions. 
References 

1. M. Traini, A. Mair, A. Zambarda and V. Vento, Nucl. Phys. A 614 (1997) 472. 

2. J. M. Cornwall, Phys. Rev. D 26, 1453 (1982). 

3. A. C. Aguilar, D. Binosi, J. Papavassiliou and J. Rodriguez-Quintero, Phys. Rev. D 
80, 085018 (2009); 

4. C. S. Fischer and R. Alkofer, Phys. Rev. D 67 (2003) 094020 [hep-ph/0301094]. 

5. D. V. Shirkov and I. L. Solovtsov, Phys. Rev. Lett. 79 (1997) 1209 [hep-ph/9704333]. 

6. A. Courtoy, S. Scopetta, V. Vento, Eur. Phys. J. A47 (2011) 49. |arXiv:1102.1599l 
[hep-ph] ] . 

7. A. Courtov. larXiv: 1107.48801 [hep-ph], 

8. A. Deur, V. Bur kert, J. -P. Chen, W. Korsch, Phys. Lett. B650 (2007) 244-248. 
[hep-ph/0509Tl3]. 

9. S. Liuti. larXiv:1108.3266l [hep-ph], 

10. E. D. Bloom and F. J. Gilman, Phys. Rev. Lett. 25 (1970) 1140; Phys. Rev. D 4 
(1971) 2901. 

11. W. Melnitchouk, R. Ent and C. Keppel, Phys. Rept. 406 (2005) 127 [hep-ph/0501217] . 

12. N. Bianchi, A. Fantoni and S. Liuti, Phys. Rev. D 69 (2004) 014505 [hep-ph/0308057| . 

13. S. Liuti, Eur. Phys. J. A 17, 397 (2003). 

14. S. Liuti, R. Ent, C. E. Keppel and I. Niculescu, Phys. Rev. Lett. 89, 162001 (2002) 

15. S. J. Brodsky and G. P. Lepage, SLAC-PUB-2447. 

16. M. R. Pennington, Rept. Prog. Phys. 46 (1983) 393. 

17. A. Courtoy, S. Liuti, in preparation. 

18. I. Niculescu et al, Phys. Rev. Lett. 85 (2000) 1186. 

19. L. W. Whitlow, E. M. Riordan, S. Dasu, S. Rock and A. Bodek, Phys. Lett. B 282 
(1992) 475. 

20. D. Amati, A. Bassetto, M. Ciafaloni, G. Marchesini and G. Veneziano, Nucl. Phys. B 
173 (1980) 429. 

21. R. G. Roberts, Cambridge, UK: Univ. Pr. (1990) 182 p. (Cambridge monographs on 
mathematical physics) 

22. R. G. Roberts, Eur. Phys. J. C 10 (1999) 697 [hep-ph/9904317] . 

23. M. R. Pennington and G. G. Ross, Phys. Lett. B 102 (1981) 167. 

24. I. Niculescu, C. Keppel, S. Liuti and G. Niculescu, Phys. Rev. D 60 (1999) 094001. 

25. R. L. Jaffe and G. G. Ross, Phys. Lett. B 93 (1980) 313. 

26. S. M. Aybat, J. C. Collins, J. -W. Qiu and T. C. Rogers, Phys. Rev. D 85 (2012) 
034043 [arXiv: 11 10.64281 [hep-ph]]. 



